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calculation of e + e — > ZZ — ► fifi + ^2/2 at and beyond LEP2 energies, where the 



We present the theoretical basis and sample Monte Carlo data for the YFS exponentiated 



left-handed part of fi is a component of an SU2L doublet, i — 1,2. The calculation is 
performed for both SM couplings and for anomalous ZZV triple gauge boson couplings 

_ _ 1 

in the conventions of Hagiwara et ai. Our formulas, which are gauge invariant, are 
illustrated in a proto-typical YFS Monte Carlo event generator YFSZZ. 
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The problem of the precision calculation of the processes e + e~ — > ZZ + 72(7) — > 
4fermions + n(j) at the higher energy range of LEP2 and at NLC (> 300 GeV) type 
energies is the subject of our paper. We are motivated to investigate this problem for 
a number of reasons, the most important of which is its role in connection with the 
verification and tests of the SU^l x U\ model of Glashow, Salam and Weinberg [Ij of the 
electroweak interaction. Indeed, at the NLC these and related processes are expected to 
play the primary role in such physics studies 0, |3|] as precision tests of the fundamental 
non-Abelian triple and quartet gauge field self-interactions in principle, for example. In 
this paper, we present and illustrate, via sample Monte Carlo data, the rigorous Yennie- 
Frautschi-Suura (YFS) [Q exponentiated Monte Carlo approach Q5| to these processes. 
The respective Monte Carlo event generator which realizes the calculation is exact in the 
infrared regime and is of leading logarithmic accuracy through 0(a 2 ) in the initial state 
QED hard radiative regime. A higher precision realization of our methods and results is 
in progress ||. 

We recall at this time the current state of the art insofar as Monte Carlo event gen- 
erators for the processes e + e~ —> ZZ + 71(7) — > Afermions + n{pf) are concerned. In 
ref. 0, both semi-analytical and Monte Carlo event generator results in the literature 
have been reviewed. What is new in our work in this paper is that we present for the 
first time the results of a calculation of theese processes in which all of the following are 
simultaneously featured: 

• realistic, finite p?, simulation of the respective multiple photon radiative effects on 
an event-by-event basis in which infrared singularities are cancelled to all orders 
in a 

• YFS exponentiated 0(a 2 ) LL initial state radiation (ISR) 
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• anomalous triple gauge boson couplings. 

Thus, our work represents the first realistic multi-photon radiative effects simulation of 
the processes under study in which the interplay between the 72(7) radiation and possible 
anomalous coupling effects can be systematically investigated. We illustrate this point in 
what follows. 

Our work is organized as follows. We first build on the recent application of the YFS 
exponentiated Monte Carlo algorithm in ref. || to the ly-pair production process by 
extending it to the Z-pair production process of interest to us here. We then present 
some sample MC data and discuss their implications for Z-pair production studies at 
high energy e + e~ colliding beam devices. Finally, we present some summary remarks. 

Specifically, in this paper, we extend the result in ref. [[UJ to the process e + + e~ — > 
Z+Z+nft) — > 4fermions+n( / y), as it is illustrated in fig. [I]. We start with the respective 
master formula. 

Referring to our master formula in ref. [|l(J , eq. (9) in this last reference, for the process 
e + + — > W + + W~ + 72(7) — > Afermions + 71(7) and to the kinematics in fig. [I], we 
arrive at the corresponding initial state YFS exponentiated 0(a 2 ) leading log (LL) master 
formula for the process e + + e~ — > Z + Z + 72(7) — > 4/ 'ermions +n{ r y) by substituting the 
respective Born level differential distribution for the latter process for the former, yielding 

1 f -A- d 3 kj f d 4 y 



da = e 2aReB + 2aBy]_ f JT ^ f 



n=0 ' j=l 



Pn(ki, . . . , k n ) 



(PpftdPpfrdSpfrdPpf. 
P Q h p\p%P% 



where we have introduced the same notation as that used in ref. |Hj for the YFS infrared 
functions B,B,D and for the YFS hard photon residuals j3 n . Thus, at the /3o level to 
which we work (i.e. f3\ and fa are included here in the LL approximation only), we 
may identify \fa with da Bom / dVLidVL 2 where dVL\(2) is the differential decay solid angle of 




Figure 1: The process e + e -> ZZ + 71(7) -> /1 + /1 + /2 + /2 + 1(7). 

/i(/2) in the respective Z's rest frame. Our respective Born cross section do Bom is taken 
from ref. [11]], where we allow for anomalous couplings gauge boson couplings f§ in the 
notation of this latter reference. We have realized the result ([!]) using MC methods of 
two of us (S.J. and B.F.L.W. |, 



12]), in complete analogy with our work in ref. [10]. 



The resulting program is called YFSZZ 1.0 and it is available from the authors [|HJ. We 
will now illustrate some of its applications. Such multiple photon radiative effects in the 
process under study here have not appeared elsewhere. 

Specifically, in fig. |2], we show the results for 10 6 simulated (weighted) events for the 
case that fx = e, f 2 = \i in fig. |T[ We show this for the NLC type CMS energy E CM = 500 
GeV. To illustrate the interplay of the n(j) radiation and the anomalous gauge boson 
couplings /4, /s, we calculate five different scenarios, /j = 0, i — 4, 5 (SM gauge boson 
couplings) at the Born level; fa — 0, i — 4, 5 (SM gauge boson couplings) at the 0(a 2 ) LL 
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Figure 2: Total cross sections and angular distributions for the process in fig. [I] for 
fx = e, f2 = // at Ecu = 500 GeV. of ot corresponds to the respective Born total cross 
section. The light solid curve, dark solid curve, solid dot, open diamond, and solid star 
correspond respectively to the differential distributions the polar angle (8) of one of the 
Z's relative to the incoming e~ direction for example for the Standard Model Born cross 
section, for the SM 0(a 2 ) LL YFS exponentiated cross section, for the fa = 0.0, /s = 0.01 
0{a 2 ) LL YFS exponentiated cross section, for the fa = 0.01, fa = 0.0 0{a 2 ) LL YFS 
exponentiated cross section, and for the fa = 0.01, fa = 0.01 0(a 2 ) LL YFS exponentiated 
cross section, respectively. 
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YFS exponentiated flo level; f± = .01, /s = 0; f& — 0, /s = .01; and, finally, f± = f$ = .01, 
where the last three scenarios are all done at the 0(a 2 ) LL YFS exponentiated /3 level. 
The respective total cross sections, along with the Born cross section denoted by of of , are 
given in the top of the figure whereas the differential cross section in the polar production 
angle of one of the Z's is plotted in the figure for each indicated. The figure 

illustrates that these anomalous couplings give angular distributions with general shapes 
similar to that of the SM couplings but that they make a significant increase in the cross 
section normalization. In this last regard, the coupling produces a stronger effect 
than does the same strength of the coupling f^, and, when both are present with the 
same strength, the effects essentially add linearly insofar as the total cross sections are 
concerned. This approximate additive linearity is consistent with the fact that f± is CP 
violating and f§ is not. At the level of the differential cross sections, there is a more 
detailed interplay between the radiative corrections and the anomalous couplings so that 
there are regions near the forward and backward directions where the effect of /s exceeds 
that of whereas in the complementary region the effect of f± exceeds that of /s, for the 
same strength of the respective couplings. Evidently, the multiple photon character of the 
radiative effects, affecting as it does the normalization and the detailed angular profile of 
the cross sections in the figure for example, must be taken into account in any precision 
analysis of such anomalous gauge boson effects. YFSZZ 1.0 allows one to do this on an 
event-by-event basis. 

In ref . |T^] , semi-analytical results are presented for the process in fig. |T] which feature 
the exact 0(a) ISR correction with soft photon exponentiation. Indeed, the authors in 



ref. [[y|] have isolated both the complete 0(a) correction and the so-called universal part 
of this correction associated with structure function evolution of the incoming e beams. 
For the process which we present in fig. [|, we may therefore compare directly the ratio the 
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(Jtot/ctot f° r the case of SM couplings, 1.081 ± 0.002, with the analogous ratio in ref. Jjj 
for the respective universal correction, which is 1.0798. The excellent agreement between 
these two numbers gives us an important cross check on our work. 

In summary, we have developed and illustrated the first ever multiple photon, finite 
Pt, amplitude based Monte Carlo event generator, YFSZZ 1.0, for the process e + e~ — > 
ZZ + 71(7) — ► fifi + / 2 / 2 + "-(7) in which infrared singularities are cancelled to all orders 
in a. We find that there is an interplay between the 71(7) radiation and the presence of 
possible anomalous couplings in the gauge boson sector. We look forward with excitement 
to further applications of our calculation in high energy e + e~ colliding beam devices. 
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Figure Captions 

Fig. 1. The process e + e" — > ZZ + 72(7) — > f x + f x + / 2 + / 2 + 71(7). 
Fig. 2. Total cross sections and angular distributions for the process in fig. [I] for f x = 
e, /2 = a ^ -^cm = 500 GeV. of ot corresponds to the respective Born total cross section. 
The light solid curve, dark solid curve, solid dot, open diamond, and solid star correspond 
respectively to the differential distributions the polar angle (6) of one of the Z's relative 
to the incoming e~ direction for example for the Standard Model Born cross section, for 
the SM <D(a 2 ) LL YFS exponentiated cross section, for the / 4 = 0.0, / 5 = 0.01 0(a 2 ) LL 
YFS exponentiated cross section, for the f± = 0.01, = 0.0 0(a 2 ) LL YFS exponentiated 
cross section, and for the / 4 = 0.01, / 5 = 0.01 0(a 2 ) LL YFS exponentiated cross section, 
respectively. 
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